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It is proved that the spectrum of scalar particles generated from the initial vacuum in inho- 
mogeneous spacetime is nearly thermal in the limit of large momentum k, if the momentum 
was defined as the variable of the Fourier transform of the coordinate in the scalar field. 
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Kandrup [ 1 ] has recently proved Parker's conjecture [2] of generating a nearly thermal 
spectrum from the background geometry of the cosmology. The proof could be straightfor- 
wardly generalized to anisotropic spacetimes. Let us give a brief sketch. 

Consider a free scalar field (p{x, t) in a spatial homogeneous spacetime, which can be 
viewed as a sum of the modes (f)k{x) Xk{t) that evolve independently one from another. 
Using a suitable time variable (say) rj, e.g., conformal time, and defining a suitable wave 
function "ifkiv) = f{v)Xk{t), where f{rj) is a function of the metric tensor, then we have the 
wave equation of "^kiv) (we leave off the index k hereafter) 



+ Q^(?7)*(?7) = 0. 



1) 



In the interval we are interested in i7 is a real, positive and finite function of the wave vectors 
k, particle mass m, curvature coupling constant ^ and metric forms g^i,. For example, in the 
case of the four-dimensional Bianchi type I spacetime with line element 



(2) 



then, after the calculation and using the conformal time variable rj = J C{t)dt, in which 
C{t) = (C^C2C3)^/^we have [3] 



^'iv) = Civ) (e ^ + + (e - ^) Riv)) + Qiv), 

^/ X _ 1 V-/ > > ^2 > _ dCi{7])/dri 
Qiri) = ^^Y^{d,-d,), d.= -^^. 

The general solution of eq. (1) can be written in the WKB form 

*(77) = (2Q)-^/2 a(r])exp{-i J''Qdri)+(3{r])exp{J\dr)) . 
For the two introduced functions a and /? we can impose one additional condition 

Q drj) — P{r]) exp{ / Q drj) 

The Wronskian condition implies that 

|a|'-|/'r = l- 



df] ^2 ^ 



(3) 
(4) 

(5) 

(6) 

(7) 



Using eqs. (5) and (6) one can reduce the wave equation (1) of ^ to two first-order coupled 
differential equations of a and /? [4,5]: 



da dil/dr) 



P exp ^2i J drf) 



(8) 
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df3 dVt/df] 



aexp (^-2i j"^ VLdrj)^ . (9) 



dr] 29. 

Note that eqs.(5), (6), (8) and (9) are exactly equivalent to eq.(l), and thus the WKB 
solution of eq. (5) is an exact form instead of an approximation. 

To make sense of the particle meaning in the initial time [rji) and final time (772), the 
spacetimes in both r]i and 772 must be either static or adiabatic static [4 ]. Suppose that the 
Fock spaces at 771 and 772 are unitarily equivalent, it then follows that ^ as /c ^ 00. 
Therefore, in the short-wavelengths limit o; — > 1 and eq. (9) becomes 

d(3 dn/dr) 



exp 



(^-2i j\d7])^ . (10) 



dr] 2Q 

If, for sufficiently large k, the ft term could be approximated as k P{ri) for a suitable function 
(3{ri), then, in a first approximation, we could integrate the above equation and obtain 

/3 = J ^^L!^exp(^-2ik J"^ bdrj^ drj. (11) 
We can define a new time variable 

T{v)^ J b{r])dr], (12) 

and then write eq. (11) as 

(3^1 F{T) exp{-2ikT) dT, 

^^^^ = wmw ^''^ 

Note that when 6(77) is a monotonic function (increasing or decreasing) then the function T 
is also a monotonic function of 77 and thus T may be regarded as a time variable. When 
6(77) is a more general function, we can cut the time variable 77 into several intervals during 
which 6(77) is a monotonic function and thus, in each interval, T may still be regarded as a 
time variable. 

For a sufficiently well behaving F{T) (this is implied from a sufficiently good behavior 
of the metric form g^„) and presuming that T ±00 for 1] —> 771^2, Kandrup then used 
the method of residues to evaluate the above integration and the thermal distribution of 
exp{— const, k) was obtained. In this case, however, it is necessary, to assume that F{t) is 
analytic in the lower-half plane except for isolated poles. 

For the case of T 7^ ±00 when r/ — > 7/1^2 we could, in the spirit of Landau's analysis of 
electrostatic damping [6], argue that for large |A;|, the integral in eq.(13) should be dominated 
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by an exponentially damped function exp{— const, k). Thus a nearly thermal form will also 
show up. 

An obviously physical interpretation for such a general character is that the short- 
wavelength modes do not probe the detailed properties of the back- ground spacetimes, 
and thus should evolve in a fairly general fashion [1]. If this interpretation is correct then 
one might expect that the thermal spectrum should also be generated from an inhomoge- 
neous background spacetime. However, the lack of homogeneity will lead the different mode 
solutions of the curved-space field equation to be mixed, and a key equation like the form 
of eq. (1) could not be obtained. Therefore, in order to study particle creation in inhomo- 
geneous spacetimes one restores the inhomogeneous perturbations of spatial homogeneous 
cosmologies [7]. As such a method could only treat the problems with small inhomogeneities, 
for studying a general situation we will in this paper directly analyze the mode-mixed equa- 
tions. It will be found that, in the limit of large momentum /c, the mode-mixed equations 
become mode-separated equations which evolve independently. Also, the general form of 
these mode-separated equations are just the form of eq. (1) and thus the thermal behavior 
arises as well. 

For convenience, we will only consider the theories in (1 + 1) dimensions, the reader could 
easily see that there is no problem to extend the models to any dimension. 

The most general differential equation describing a scalar field in an inhomogeneous 
spacetime is 

f + Ei., § + Z,(., «) |i + C(., t) g + B(., t) ^ + Ai., ^ 0. (14) 

where A{x,t), ...E{x,t) are functions of the particle mass, curvature coupling constant and 
metric forms, which need not to be specified in our proof. 

Defining Ak{t), Ek{t), and (j)k{t), as the Fourier transforms, with respect to the variable 
X, of A{x, t), E{x, t) and 0(x, t) we can obtain from the above equation 

<l>k- {k-lfEi<f>k-i-i {k-l){Di(f>k_i+CicPk-i)dl+ {Bicf>^_i+Ai(f>k-i)dl^O, (15) 

J—oo J—oo J—oo 

where an overdot means d/dt. This differential-integral equation manifests mixing between 
the different modes as a consequence of the inhomogeneity of the spacetimes. 

Using the discussions described below eq. (13) we know that the exponential damping 
factor exp{— const. I) should dominate in the limit of large / for the functions Ai{t), Ei{t),. 
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Thus, in the first approximation of large A;, we can reduce the above equation to 

'^,-k^E^k-tk{Di, + C^k) + {Bik + A^k) =0, (16) 

where 

£;= / Eidl, D= Didl, C= Cidl, B= Bidl, A= Aidl, (17) 

J —oo J —oo J —oo J —oo J —oo 

One then finds that in the hmit of large momentum k the mode-mixed evolution equation 
becomes a mode-separated evolution equation in any inhomogeneous spacetime. 

Furthermore, as k is very large, we can, in the case of sufficiently well behaving functions 
A{t), ...,E{t), approximate eq. (16) by 

4>k-k^E4>k-ikD4>i^ =0, (18) 

After defining the wave function "if kit) = fit) <i>k , in which 

hit) = exp (^-\ikj D{t) dt^ , (19) 

we then obtain the wave equation of "if kit) 

tfe(i) + ^(t) = 0, = k^ - E^ . (20) 

As the above equation has the same form as the key equation, eq. (1), we could thus use 
the same procedure as described in the earlier section to prove that the particle spectrum 
generated from initial vacuum in inhomogeneous spacetime is nearly thermal in the limit of 
large momentum k, if the momentum was defined as the variable of the Fourier transform 
of the coordinate in the scalar field. 

Note that one needs not worry about the situation of raising a non-positive value of 
Q^(i) in the last equation. In fact, only for the cases with reasonable metric forms and after 
choosing a suitable time variable could one make sense of the particle meaning in the initial 
and final times, which are necessary conditions for the study of particle creation in curved 
space. Under such a condition one sees that the particle will propagate with a wave form, 
and thus the value of the f2^(t) term in the last equation must always be positive. 

Finally, we want to remark that the scale above which the spectrum is exponential may 
be determined by the characteristic length over which the universe is inhomogeneous. The 
spectrum can at best be thermal only for wavelengths short compared to this characteristic 
length. 
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